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The use of a running coupling constant in renormalizable theories is well known, but the imple- 
mentation of this idea for effective field theories with a dimensional coupling constant is in general 
less useful. Nevertheless there are multiple attempts to define running couplings including the effects 
of gravity, with varying conclusions. We sort through many of the issues involved, most particularly 
the idea of operator mixing and also the kinematics of crossing, using calculations in Yukawa and 
\(f) 4 theory as illustrative examples. We remain in the perturbative regime. In some theories with 
a high permutation symmetry, such as \4> 4 , a reasonable running coupling can be defined. However 
in most cases, such as Yukawa and gauge theories, a running coupling fails to correctly account for 
the energy dependence of the interaction strength. As a byproduct we also contrast on-shell and 
off-shell renormalizaton schemes and show that operators which are normally discarded, such as 
those that vanish by the equations of motion, are required for off-shell renormalization of effective 
field theories. Our results suggest that the inclusion of gravity in the running of couplings is not 
useful or universal in the description of physical processes. 



1. INTRODUCTION 

Quantum corrections to scattering processes include a 
kinematic dependence on the energy scale of the process. 
In renormalizable theories, the idea of a running coupling 
constant absorbs a dominant and universal set of quan- 
tum corrections into a well defined logarithmic function 
of the energy, making this the appropriate expansion pa- 
rameter in perturbation theory. Physical processes at 
a given energy scale are best expressed in terms of the 
running coupling constant defined at that scale. 

Quantum corrections due to the gravitational interac- 
tion to various processes are also calculable using effective 
field theory methods [lj. Because the gravitational cou- 
pling carries a dimension, the quantum corrections to a 
matrix clement Ai carry a power-law dependence on the 
energy scale 

M~a[g + b 9 K 2 q 2 + cgn 2 q 2 log(-q 2 ) + ...] , (1) 

where g generically denotes a coupling constant or com- 
binations of constants and 

K 2 = 327rGNcwton = -TTJ > (2) 
IVl p 

with G being Newton's constant and Mp = the Planck 
mass. Following the success of running couplings in other 
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contexts, it is tempting to try to also absorb some of the 
gravitational corrections into a running coupling constant 
g(q 2 )- We note many such attempts [3|-|15|. 

It is by now well known that application of the renor- 
malization group to effective field theories such as grav- 
ity does not lead to a traditional running coupling con- 
stant. In his influential paper on effective field theory 
, Weinberg showed that the content of the renormal- 
ization group in these theories is to relate the highest 
powers of q 2 log(— q 2 /fi 2 ) to each other. This behavior 
has been explored subsequently in more detail [l7rl2l| . 
This is due to the power counting relations of effective 
field theory which tell us that loop processes generate 
higher order operators that involve more powers of the 
derivatives and/or fields. Because of the increasing pow- 
ers of q 2 in the factors of q 2n log"(— q 2 //i 2 ) this does not 
lead to a renormalization of the leading coupling con- 
stant, but rather the higher order couplings are renor- 
malized. In addition, the logarithms enter differently in 
different processes or even in two form-factors for the 
same process^ 1 . 

Therefore attempts to define a running coupling nec- 
essarily involve definitions which fall outside of the usual 
renormalization group. They tend to involve attempts to 



Because the coefficients of the higher order operators are the 
ones that absorb the divergences of the effective theory, they do 
have a dependence on the scale fi that occurs in dimensional 
regularization, but this dependence does not induce a running in 
the original lower order couplings. 



2 



identify a typical q 2 in an amplitude with a renormaliza- 
tion scale 2 M 2 and absorbing the higher order q 2 effects 
into the coupling constant 

g(M 2 ) =g + bgn 2 M 2 + cgn 2 M 2 log(M 2 ) . (3) 

It is worth exploring whether such a definition can make 
sense in physical processes. 

This procedure is also outside of the usual effective field 
theory methods. One applies effective field theory as an 
low energy expansion about zero energy. Higher order 
energy dependence is associated with operators that are 
higher order in the energy expansion 3 and does not renor- 
malize the low order operators of the theory. A coupling 
constant definition such as Eq. [3] attempts to define the 
theory around the high energy renormalization scale M. 
In this case, success could be achieved if the definition 
correctly captures the effects of the loop corrections at 
this energy scale. 

Some potential pitfalls are visible in this strategy. 

• One is that in the effective field theory context, 
higher momentum dependence is associated with 
new operators carrying extra derivatives, not the 
original operator carrying the coupling g. This 
brings in the need for operator mixing, as renor- 
malization conditions at a given scale can involve a 
mixture of a set of operators. 

• However, even if one can define a running coupling 
using some combination of the relevant operators, 
there is the question of universality. That is, one 
has to assess whether the operator mixing is the 
same for all processes. It is possible that a defini- 
tion that is useful for one reaction may be deleteri- 
ous when used in another reaction. 

• In addition, there is a problem of kinematics and 
crossing, in that the energy variable q 2 can take on 
different values and even different signs in different 
contexts. For example, in space-like versus time- 
like reactions, q 2 changes sign, so that a running 
coupling that decreases with energy for a space-like 
process will increase with energy if that same reac- 
tion is crossed into the time-like regime. Crossing 
will turn out to be a major obstacle in the Yukawa 
and gauge theory cases. 

In this paper we will explore these issues using calcula- 
tions of the gravitational corrections to A</> 4 and Yukawa 
theories. We use these theories as test cases in order to 



2 In this paper we will use the notation E or M for various defini- 
tions of the renormalization scale and reserve /i for the scale that 
arises in dimensional regularization through the factor /i 4_d in 
front of dimensionally regularized loop integrals. 

3 If particle masses are considered, there is some renormalization 
of low order operators, however, we consider massless theories in 
this paper. 



avoid the irrelevant complications of gauge invariance. 
However the general lessons of our results will apply to 
other theories also. We will find that outside of some 
special cases, the idea of the gravitational contribution 
to the running of a coupling constant is not a useful idea 
in the perturbative regime. 

2. PREVIEW OF KEY ISSUES 

Let us first review various ways of calculating the run- 
ning coupling in renormalizable theories. 

At the most physical level, one can calculate any given 
physical process including quantum corrections and iden- 
tify the large logarithms that can be absorbed into the 
running coupling. However, the use of physical processes 
can sometimes be complicated by the presence of imagi- 
nary parts to the amplitudes and by not being sure what 
part of the quantum corrections are universal enough to 
be absorbed in the coupling constant. 

To address these issues, one can alternatively define 
the coupling by renormalizing at an unphysical Euclidean 
point p 2 — —M 2 , avoiding the cuts and poles of the phys- 
ical amplitudes. 

A third method of great practical utility is to study the 
divergences of the coupling constant. In dimensional reg- 
ularization, the I/e divergences are always accompanied 
by log^i, where a factor of // is introduced to keep the 
dimensionality of loop integrals unchanged. By dimen- 
sional analysis then, if the only large scale in the theory 
is the renormalization scale M, the log fj, dependence al- 
ways tracks the logM dependence. The 1/e behavior is 
obviously universal, since it goes into the renomalization 
of the coupling constant, and the accompanying logM is 
also universal and can be readily incorporated into the 
running coupling constant g(M). This is powerfully ex- 
ploited in renormalization group arguments to show that 
this running coupling is the appropriate coupling for all 
processes at this scale. 

If we look at the nature of the gravitational correc- 
tions, we see some crucial differences. The well-known 
presence of divergences in gravitational loops is not itself 
a significant issue, but the nature of the renormalization 
procedure dealing with them is important. Because the 
gravitational coupling carries inverse mass dimensions, 
the divergences go into the renormalization of new oper- 
ators that carry extra derivatives. For example, when we 
discuss the Yukawa couplings of a scalar 4> and fermion 
ip, originally of the form 

C Y = -gW>ij) ■ (4) 

At low energy, the divergences go into the renormaliza- 
tion of the coefficients of the higher order operator 

O g3 =g 3 ^d^. (5) 

or, as we discuss below, into higher dimension four- 
fermion operators. Following the 1/e behavior of the 
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loops will not capture the renormalization of the origi- 
nal operator. 

Even if we give up this useful technique, we can still 
study the energy dependence of physical processes or con- 
sider the strength of the interaction at a Euclidean point 
p 2 = —A/ 2 . It is clear that this will then involve a linear 
combination of the initial operator plus the higher order 
operators. This is what we mean by operator mixing. 
A rule for such a procedure can always be developed to 
define a coupling at a given energy and this will yield a 
given definition of a running coupling that includes power 
law running. But the question then is whether this defi- 
nition is useful. To be useful, it should be in some sense 
universal, so that it applies to other reactions also, and it 
should encapsulate at least some of the large corrections 
to physics processes. The presence of multiple operators 
in the effective field theory basis argues against universal- 
ity - different operators contribute differently from pro- 
cess to process. And power-law kinematics also argues 
against utility, because as mentioned previously a kine- 
matic variable that is positive in one reaction is negative 
in a related reaction. A coupling that minimizes the en- 
ergy dependence of one process with increase the energy 
dependence of the related process. 4 

It is possible that this procedure can still work. In 
A</> 4 theory we will see that physical processes regularly 
involve a mixture of spacelike and timelike subdiagrams, 
because of the permutation symmetry of the original in- 
teraction. In this case, we will be able to define a reason- 
able running coupling with power law running, subject 
to only modest ambiguities due to the renormalization 
scheme. 

However, in most theories we find that the power law 
running is not a useful concept. In the perturbative 
regime, these theories are better described by an oper- 
ator basis with coupling constants that do not run nor 
mix. 



GRAVITATIONAL CORRECTIONS TO A0 4 
INTERACTION 



In this section we will explore the various ways of defin- 
ing a running coupling in \(f) theory. This effort is rea- 
sonably successful, and provide an illustration of what 
powerlaw running could look like. The feature that is 
most important in this construction is the mix of space- 
like and timelike diagrams, with a high permutation sym- 
metry. 

We consider a massless real scalar <fi with a A0 4 in- 
teraction, coupled minimally to gravity. The Lagrangian 



reads 



A ,4 

4!* 



(6) 



where k 2 = 327rGNowton, 9^ is the metric tensor, R is 
Ricci scalar, and A is the scalar-self coupling. 

Temporarily ignoring the gravitational interaction, the 
one loop scattering amplitude in this theory is derived 
from the diagrams of Fig. [T](a), (b), and has the form 



iM = + 



3U2_ 
32^ 

lo 



h log 47T — 7 



2 (4irY 
~ l0g (? 



(7) 



where e = 4 — d and A is the bare coupling constant, and 
the channels s, t and u are defined as usual s = (P1+P2) 2 , 
t — (Pi ~ Ps) 2 : u — {P2 — Pz) 2 - In order to use the on-shell 
process to define a running coupling, we can choose to 
measure the renormalized coupling at the point s = 2E 2 , 
t = u = —E 2 . This lets us define the effective coupling 
constant \{E) as 



i\[E) 



-iX + 



3i\ 
32^2 



2 r 



log Air — 7 



2(4tt) 



log 



2E 2 
V 2 



2 log 



M 2 



(8) 



such that the on-shell perturbative scattering amplitude 
becomes 



-iM = -i\(E) - 



i\ 2 {E) 
2 (4tt) 2 

+ lo 



log( 



2E 2 ) 



log 



— u 
E 2 



E 2 ) 
(9) 



For all s, t, u of order E 2 all the logarithms are small. 
The potentially large logs have been absorbed into X(E). 
The quantum corrections proportional to A 2 vanish at 
the renormalization point and are small throughout the 
physical region. The j3 function is calculated from 



B(\) E dX{E) 



3A 2 
16tt 2 



(10) 



It is often common to rcnormalize a symmetric off- 
shell Euclidean point. In this case we treat all lines as 
incoming and choose kinematics p 2 — —M 2 , s = t = 
u = —AM 2 /i. This allows a definition 



4 Note that the sign of the momentum is not an issue with a run- 
ning coupling with logarithmic behavior, log(— |q 2 |) = log(|q 2 |) + 
i-K. The magnitude of the logarithm is universally present in both 
spacelike and timelike processes, while the imaginary part is part 
of the residual quantum correction. 
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(11) 
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with a similar expansion of the matrix element, and re- 
lated beta function. Finally, we note how the beta func- 
tion can be read off from the coefficient of the 1/e, since 
it is intrinsically connected to the log fj, in the scattering 
amplitude. 

We now turn to the comparable definitions of the run- 
ning coupling with the inclusion of gravitational correc- 
tions. We employ two different methods. In all methods, 
we use dimensional regularization scheme to regularize 
our integrals. First, we define the running coupling con- 
stant as an effective coupling for scattering processes. In 
the second method we use an off-shell procedure to cal- 
culate the /3 function. These will yield similar results. 

3.1. Gravitational corrections in on-shell scattering 
processes 



In principle, one needs to include as well the scalar masses 
and the contribution from soft gravitons that potentially 
will remove the IR divergence. We do not follow this pro- 
cedure here. Instead, we choose to ignore altogether the 
contributions from C (x) . Another procedure that avoids 
these divergences will be followed in the next sections, 
where we perform all the calculations off-shell. 




P2 / \Pt 



+t + u ' +3 others 

(a) (b) (c) 



Let us now include gravity. To find the graviton propa- 
gator one perturbs the metric tensor about the flat back- 
ground g^y = T}^ LV + k/i m „, where h^ v are the spacetime 
fluctuations. Then, one expands R in terms of and 
writes R in the form h^ v 0^ v,f " T h pa . To simplify the cal- 
culations, we fix the gauge freedom by employing the 
harmonic gauge d^h^ — d u h"/2 = 0. Finally, one ob- 
tains the graviton propagator 



M 2 



(12) 



We will first explore an on-shell renormalization 
scheme. We consider the different diagrams contribut- 
ing to the on-shell scattering process <fi + cf) —> 4> + 4>. In 
addition to the quantum scalar corrections to the s, t and 
u channels (Fig. Q](b)), we include the gravitational cor- 
rections to the wave- function and vertex, as shown in Fig. 
Q] (c), (d), (e) and (f). The gravitational wave-function 
renormalization and diagram (e) vanish for massless par- 
ticles. Also, diagram (d) does not contribute to the ver- 
tex corrections in the dimensional regularization scheme. 
Hence, we are left only with diagram (f) 



A 



U) 



IK X 

-u lo 



slog 



<log 



j [s 2 C(s) + t 2 C(t) + u 2 C(u)] , (13) 



2(4tt) 

where fj, is introduced in the scalar vertex as /i 4-d A<^> 4 /4! 
for the purpose of dimensional regularization. The func- 
tion C{x) is given by the integral 



C(x) 



dwdt 1 . (14) 
o Jo «>£(! ~0 X 




(d) 



+3 others 
(e) 



+5 others 

(f) 



FIG. 1: The contributing diagrams to the running of the 
scalar coupling constant A. 

Note that there is no UV divergence in the above on- 
shell scattering amplitude. This is a special feature of 
(j) 4 theory and arises from the permutation symmetry of 
the Feynman diagrams. While individual diagrams are 
of course divergent, the divergent pieces sum to 



Mdiv ~ Xk 2 -(s + t + u) . 

e 

However this sums to zero since we have 



t + u = 



(15) 



(16) 



when evaluated on-shell 5 . Note also the related feature 
that despite the apparent presence of fi in the scattering 
amplitude Eq. 1131 in fact this amplitude is independent 
of fx because of the on-shell identity Eq. [THJ 

A related unique feature for this theory is that the 
higher order operator vanishes by the equations of mo- 
tion, or equivalently can be removed by a field redefini- 
tion. The higher order operator that is generated by one 
loop gravitational corrections is 



(17) 



This would generate a matrix element proportional to 
s + 1 + u, which as we have seen vanishes. By integrating 



These integrals are IR divergent which makes their inclu- 
sion problematic in calculating the scattering amplitudes. 



J If we had used massive scalars, we would obtain UV poles ~ 
m 2 (l/e — log(p 2 //^ 2 ) + finite). 
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it by parts it can be seen to be equivalent to the operator 
</f 3 D^>, which vanishes by the equation of motion. Such 
operators can be removed from the operator basis in favor 
of other local operators. In this case we can see this by 
performing the field redefinition 



Ai 
3 ' 



(18) 



which removes the term C\ t and to first order in Ai gen- 
erates the operator ~ XX\4> 6 , which comes from the ex- 
pansion of the original tfi 4 interaction. In dimensional 
regularization this cf> 6 operator does not mix with tp 4 at 
one loop since the massless tadpole loop integral van- 
ishes. Therefore for onshell renormalization we do not 
need to consider operator mixing at this order. We will 
however return to this issue in the next section. 

Now, we choose to renormalizc (fT"3|) at the physical 
point s — 2E 2 , t = u = —E 2 , and then define the effec- 
tive coupling constant X(E) as 



i\(E) = -iX - 



3iX : 



2(4tt) 



ibgOVA 2 

9 Cj i 

(4^) 2 



(19) 

where the second term is the quantum scalar correction, 
and we have omitted an imaginary phase. Now the full 
amplitude is given by 
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~ 2 (4tt) 2 

+ulog 



—u 
'E 2 



in 2 X(E) 
' 2 (4tt) 2 



s 2 C(s) +t 2 C{t) +u 2 C(u)] , (20) 



This has been a useful definition as the hard quantum 
corrections for the gravitation loops vanish at the renor- 
malization point and stay small in the physical region. 
Using this definition, the /3 function reads 



log(2)« 2 A f2 
8n 2 E 



(21) 



Under this procedure, the gravitational corrections do 
not tend towards an asymptotically free theory while in 
the perturbative region. 



3.2. Off-shell renormalization 



no essential complication, because we know all the oper- 
ators that can be renormalized. It actually provides a 
simplification, as we can choose a convenient symmetric 
point, and in the Euclidean region can avoid poles and 
cuts. However, in our case, the off-shell point brings in 
a potential complication, as new divergences appear and 
the higher order operator of Eq. [T7] now gives a non- 
vanishing matrix element. These two issues are related, 
as the divergence is absorbed in the coefficient of the 
higher order operator. Thus we have operator mixing 
appearing at the off-shell point, while it did not appear 
on-shell. However, by dealing with this feature we can 
still obtain results similar to the on-shell running cou- 
pling. 

In this method we compute the scattering amplitude 
using off-shell momenta. In this case, the diagrams of 
Fig. [T] generate the operator 0\ 1 = — Ai0 2 9 A1 09 M </)/8. 
The scattering amplitude is given by 
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(22) 



and, for example, 

C(pIpI,s) 
i pi 



dwd^ 



o (i-«0(p?(i-O+p1O + <(i-O« 



(23) 



In this section, we provide another method to calcu- 
late the (8 function by using an off-shell renormalization 
point. In renormalizable theories, going off-shell provides 



We then compute the scattering amplitude at the Eu- 
„2 _ „2 _ p 2 _ p 2 _ _^fi_ and 



clidean momenta p\ 



P\ 



6 



-4M 2 /3. Hence, we find 



3.3. Lessons from cf> theory 
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where p = — i(Air) 2 Z / n 2 \M 2 is some numerical coeffi- 
cient whose value does not affect the f3 function 



1 log [4^(1 -Q/3] 
* -3 + 4£-4£ 2 



2.83. 



(25) 



In order to deal with the issue of operator mixing, we 
need to choose appropriate renormalization conditions. 
Now, we define effective couplings X(M), and Ai(M) such 
that the scattering amplitude is given by 



A = —iX(M) + iM 2 \i(M) 
For the higher order operator we define 
dA 



dM 



= iX 1 (M). 



(26) 



(27) 



Because Ai is generated through loops and is higher order 
in the energy expansion, we treat it as a quantity of order 
A 2 or Xk 2 . This allows us then to neglect any feedback 
from Ai into the original coupling A. Such feedback would 
occur from further loops and would be of higher order 
AiA ~ C(A 3 ) or Xik 2 ~ 0{X 2 k 2 ), both of which we drop. 

Using these definitions, the divergence ~ 1/e is ab- 
sorbed as usual in A, and the new divergence ~ (s + t + 
u) /e is absorbed into Ai . We then solve the system (|26p 
and (|27p simultaneously to find 



A(M 2 ) = A + 
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log 



2 (Any 
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(28) 



Finally, the f3 functions of A and Ai read 
r dX(M) 3A 2 k 2 A 



/3(A) 
/3(Ai) 



71/ 
71/ 



M 2 

dM ( 47r ) 2 4tt 2 
d\i(M) _ 3A 2 k 2 X 
dM ~ (4tt) 2 M 2 4^" 



(29) 



We have used two different methods to define the run- 
ning coupling. Both on-shell and off-shell Euclidean 
methods yielded similar results. The higher order opera- 
tor in this situation vanishes for on-shell matrix elements, 
so that the on-shell method did not require any opera- 
tor mixing and resulted in a finite amplitude. However 
we needed this unphysical operator and operator mixing 
in order to accomplish the renormalization when work- 
ing off-shell. Of course when we continue and apply this 
operator to a physical process it will again vanish. The 
slight differences in the beta functions can be accommo- 
dated by a scheme dependent renormalization scale. 

In renormalizable theories, we also know that the run- 
ning coupling will also be the one relevant for loop di- 
agrams at a given energy. One way to show this is 
to use renormalization group arguments. However, we 
also know this from the structure of Feynman diagrams. 
While loop momenta run over all energies, in dimensional 
regularization the loop integrals are dominated by the 
overall energy scale of the problem, aside from infrared 
and collinear regions that can be dealt with using other 
means. This feature also makes sense given the progress 
in constructing loop results from unitarity cuts, which 
use the physical on-shell amplitudes. 

We have not used this running coupling in a calcula- 
tion involving higher orders in the loop expansion, al- 
though we expect that this coupling remains an accept- 
able one. As mentioned in the introduction, in this case 
the renormalization group does not dictate the utility of 
the running coupling. However it appears likely that this 
coupling will appear in higher order processes. Because 
of the simplicity of this theory, the <fi 4 interaction is the 
only one involved in higher order interactions. Within 
dimensionally regularization of loops for processes at an 
energy E, the only relevant momentum scale is again is 
the energy E because the particles are massless. Disper- 
sive techniques will use the on-shell amplitude, and hence 
will involve the coupling that we defined initially. Wick 
rotation of Feynman amplitudes would transform ampli- 
tudes to Euclidean momenta, where we found a similar 
result. 

We conclude that this definition of a running coupling 
is a useful one in tp 4 theory at one loop order and may also 
be useful at higher order in perturbative calculations. 



4. GRAVITATIONAL CORRECTIONS TO 
YUKAWA INTERACTIONS 

In this section we follow the same lines above to calcu- 
late the j3 function of Yukawa interactions. Here the ba- 
sic vertex cannot be defined with all legs on-shell. Much 
as in gauge theory, on-shell renormalization requires a 
scattering amplitude with two vertices (the equivalent of 
Coulomb scattering in QED). However the off-shell func- 
tion vertex can be defined at an unphyiscal kinematic 
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point. We again explore both on-shell and off-shell renor- 
malization. 

If we define the original Yukawa coupling constant by 
the Lagrangian 

C Y = T<jnpi> , (30) 

then we are looking for a running coupling of the form 

T(M) = T + aTK 2 M 2 (31) 

for some constant a, when implemented with a renormal- 
ization scheme defined at the scale M. It will be possible 
to make such a definition. 

However, we will also find that any such definition 
does not correctly capture the loop effects of the quan- 
tum corrections in all relevant processes. Let us highlight 
the main issue here before providing the explicit demon- 
stration. Consider two physical processes involving the 
Yukawa couplings such as // — » // and // — > // which 
proceed through the exchange of the scalar field. These 
of course are related by crossing, with the momentum- 
squared, q 2 , of scalar changing from positive for the 
time-like process to negative for the space-like process. 
When calculated explicitly, the loop effects from the ver- 
tex function in these scattering processes will depend on 
the variable q 2 . In this case the loop corrected matrix 
element will have the form 



M 



r(MWr(M) + loops 



(r 2 + 2aT 2 K 2 M 2 + 2a'T 2 K 2 q 2 



, (32) 



where a' is a number that emerges from the loop calcu- 
lation. The problem is that even if the definition of the 
running coupling is chosen in such a way as to capture the 
main quantum corrections for one process, say // — > //, 
it will have the opposite effect in the crossed process. No 
definition of a running coupling can summarize the quan- 
tum corrections in both processes because the quantum 
effects go in different directions in the two cases. If they 
make the matrix element smaller in one channel, which 
naively looks like asymptotic freedom, they make the am- 
plitude larger in the other process, which does not look 
like asymptotic freedom for the coupling. 

On the other hand, operator mixing with a higher di- 
mension operator does correctly describe the quantum 
effects in both channels. Because the factors of q 2 cancel 
in the loop effects 



(2aTVg 2 )l - 2aTV 



(33) 



these effects are described by a contact operator 

a'T 2 K 2 ipiprpip ■ (34) 



4.1. Operator mixing 

Gravitational corrections to the vertex (as shown in 
Fig. [5]) will generally generate the higher order operators. 
A convenient basis for our calculation can be chosen to 
be 



2 = cb (tpd 2 ^ + d 2 ^) , 

3 = 4>d^d^. 



(35) 



and include their respective coupling constants 

C ho ^g 1 1 +g 2 2 +g 3 3 . (36) 

Despite appearances, all three of these can be shown to 
vanish by the equation of motion. For 0\ this follows 
from relating the operator to 



4 = <i>d^ri v d u ^ , 



(37) 



which clearly vanishes by the Dirac equation, and then 
using the identity 



(38) 



which follows from expressing the LHS in terms of com- 
mutators and anticommutators. This turns O4 into a 
combination of other operators 



04 = 03 - iOi 



(39) 



In addition O3 can be seen to vanish through equations 
of motion through integration by parts. Here we define 



we find that 



5 = {d 2 <f) H . 



05 = 20 3 + 2 



(40) 



(41) 



However, a little more care is needed in this discussion 
because all legs of this three-point vertex cannot be on- 
shell at the same time. It is easy to see that an operator 
such as 03 can lead to a non- vanishing matrix element 
in a physical process, such as // — > //, where the scalar 
field is off-shell. However, this matrix element is a con- 
stant independent of momentum, and is equivalent to a 
local four-fermion operator. What this means is that an- 
other set of operators needs to be introduced which are 
used for onshell processes instead of the operators 0j 
which vanish by the equations of motion. These are 

Qi — xpip^tl) 

Q 2 = ^ v ^g^ 

Q 3 = 4>8,Md^ 

(42) 



This works for both processes as the answer is indepen- 
dent of the sign of q 2 . 



Direct calculation using Oi in tree level physical processes 
shows that their effect are equivalent to contact operators 



such as the Qi listed above. Equivalently one can use field 
redefinitions to remove the Oi operators, and these can 
generate the contact operators. 

If we define the couplings of operators of Eq. HUby 
we will sec that the qi will be infinitely renormalized by 
loop processes, and hence they are needed in the descrip- 
tion of the on-shell renormalization. 



(a) 



(b) 







4.2. Off-shell renormalization 

We first consider the vertex function directly and the 
quantum corrections to it as shown in Fig [5J Because we 
need to treat this vertex off-shell, we will need to include 
the unphysical operators that vanish by the equations of 
motion Oi. We will see that it is possible to extract the 
effects of each of the three operators, using a set of renor- 
malization conditions. Our treatment is relatively brief 
as the methods are similar to the off-shell renormalization 
of tfi 4 theory. 




FIG. 2: 

theory. 



Gravitational corrections to the vertex in Yukawa 



In this method we use the Euclidean momenta p\ = 
— £iM 2 , and p\ ■ pi = —&M 2 , where p\ and P2 
are the fermions momenta, and £1 and £2 are arbitrary 
positive parameters. Unlike the case of A<^ 4 theory, where 
we average over the different scattering channels, the loop 
results in the case of Yukawa theory will generally depend 
on the kinematics of the problem. In the following we 
parameterize this dependence using the parameter 77. For 
time-like energy variable q 2 we have rj = 1, while for 
space- like q 2 we have 77 = — 1, where q = p\ + r\p2- 

In addition to the vertex diagrams shown in Fig. [2] 
we need to calculate the self-energy diagrams of Fig. 
[5J The self energies and vertex corrections acquire log- 
arithmic dependence of the form log (£iM 2 //i 2 ) , and 
log + £,i)M 2 / p?) . We immediately see that only 
£1 = 0, or £1 = —77^2 are problematic, and hence we 
avoid these values in the following analysis. We find that 
the total scattering amplitude, apart from finite pieces 
that do not affect the /3 functions, is given by 



FIG. 3: The diagrams contributing to the self energies of (a) 
fermions, and (b) bosons. 



A = -ig + iM 2 (2£i52 -653)+ giPi^o^ 



Si 



S-> 



8 (4tt) 2 [ £1 - rfa 

4^i-^4(a + ^ 2 )M 2 s 2 

2 (4tt) 2 4 (4tt) 2 



P\^P2^ 



4 (4tt) : 



(6+2r ? 6)M 2 l S 1 , 



(43) 



where 
Si 
S 2 



2 , „ f^ M2 

7 + log 47T - log 7T— 

e V /» 

- - 7 + log 4tt - log 

e V M 



(44) 



In order to accomplish the renormalization we will de- 
fine a set of renormalization conditions. We use the kine- 
matic and Dirac structure to isolate the relevant terms. 
The overall amplitude is defined as 



A 



-ig(M) + 1M 2 (2&fla(M) - MM)) 
- 9l {M) Plll p 2v a^ . (45) 



By inspection, we find, apart from a trivial additive con- 
stant, 



9l{M)=g x 



V9 K , (M* , 
T2 lo S ( —T ) . ( 46 ) 



8 (4tt) z 



where we have absorbed the pole ~ Px^Piv^ 1, '/e in <7i- 
We can also define g-z{M) and gz{M) through the defini- 
tions 



dA 



d&M 2 ) 
dA 



2ig 2 (M) 



(47) 



5(&M 2 ) 

Hence, solving (|4"5"|) and (|4T|) we find, apart from trivial 



9 



additive constants, 
9(M) = g 



5g 3 f^M 2 
2 lo § 



2(4tt 



4 (4tt) 2 



92 (M) 
53 (M) 



5 

52 - t 



5 



4(47r)^!M 2 



53 



gn 2 r] 



4 H ^ 



-log 



M 2 



(48) 



where the pole ~ 1/e is absorbed in gr, the pole ~ M 2 /e 
is absorbed in 173, while g% does not get any pole contri- 
bution. 6 Finally, the (3 functions read 



P{92) 

/3(53) 



5g 3 



(4tt) 2 2 (4tt) 2 
rjgK 2 



4(4tt) z 

55 3 
(4tt) 2 CiM 2 

TjJK 2 



2 (4tt)' 



(49) 



Hence, we see the gravitational correction of [3(g) de- 
pends on the kinematics through the parameter r\. More 
on this is discussed below. 





FIG. 4: Tree diagram for the on-shell scattering processes 
involving fermion.The filled circle denotes the set of vertex 
renormalization diagrams. 



former has only a t-channel exchange and the latter only 
s-channel. Because these processes are on-shell, we can 
drop the operators 0; associated with the vertex itself, 
but must include the four fermion operators Q\ and Q2 
associated with the four-fermion process. 

The key diagrams occur via the exchange of a scalar 
boson and so include the vertex correction on either side 
of the diagram, as in Fig. @] This set is the analogous to 
the set of diagrams considered for the running coupling 
in renormalizable theories. The process also includes a 
set of other diagrams, shown in Fig. [5] While we have 
calculated the divergences in these diagrams and veri- 
fied that they can be absorbed in the coefficients of the 
four fermion operators Qi, we do not include them in the 
definition of the running coupling. 




+3 others +3 others 

FIG. 5: Box diagrams. 

For this calculation we can use our previous results 
with the on-shell condition £1 = 0. We find that the 
gravitational corrections to the self energies vanish as it 
goes as ~ K 2 pf = 0. For the on-shell process // — > // 
we find that the matrix element including gravitational 
corrections is given by 



iM = A, 



on shell— a " 



(pi + vp 2 ) 



on shell— & j 



(50) 



where 



A 



on shell— a 



-1-9 a 
1 



^75a K 
(47T) 2 

3 



2 Sl 16 



P1-P2- 



(51) 



4.3. On-shell renormalization 



and 



Finally we turn to the renormalization of the Yukawa 
coupling through an on-shell process. Similarly to the on- 
shell renormalization of the gauge couplings in QED, we 
consider the one-loop corrections to a scattering process 
such as / + / -> / + / or / + / -> / + /. For clarity 
in separating the crossed channels, we will in this section 
refer to two flavors of fcrmions f a and ft, so that we will 
compare f a + f b -> f a + f b or f a + f a -> f b + f b where the 



6 This can be seen by taking the derivative of (143 1) with respect to 
Z1M 2 . 



Si 

S2 (Pi -P2) 



7 + log 47T — log 

e 



7 + log 47T — log 




(52) 



In the above analysis we have restored the fermion mass 
in order to avoid an IR singularity in the logarithm. No- 
tice that the result ([ST]) can be obtained from (l43l) by 
setting £i=0 everywhere except inside the logarithm 
where it is replaced by the fermion mass, and replacing 
M 2 —> —M 2 , i.e. working with physical rather than Eu- 
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clidean momenta. Hence, 
-iM(pi -p 2 =&M 2 ) 

\ ig a gbK 2 



{ 2 (47iT 
+ (ao b)} 



I9a9b 

' ~ 2r/£ 2 M 2 



3_ 

Tg 



(53) 

and we have used pi ■ p 2 = £,2^1 2 . 

In order to use this result for the renormalization of g, 
we define the i-channel process f a + fb — > f a + fb with 
our kinematic features defined by r\ = — 1, £2 = 1/2, 
q 2 = —2pi ■ P2 = —M 2 defining the matrix element 

- iM{f a + h^fa + fb) = HsD-o Hs&) 

q 

+ (a b) - iqi . (54) 

Here we have included only the coefficient of the four- 
fermion operator Qi = i/} a ip a ipb4>b, as described in Eq. 
1421 The other operators are not needed for this analysis. 

We can absorb the pole in (1531) ~ 1/e, which is due to 
gravitational corrections, into qi by defining 



<7i = <7i 



g a gbn 
(An) 2 



M 2 



1 



(55) 



The residual loop effects are contained in the renormal 
ized value of g and we find 

3g l n 2 M 2 



9i 



(56) 



32(4tt) 2 
for i = a,b. 

At this point, one can calculate the gravitational con- 
tributions to the P function of g and q\ as 

3gtK 2 M 2 



P(9i) = M 



M 

dM 



and 



Hyq J dM 



16(4tt) 2 

3g a g b K 2 
4(4tt) 2 



(57) 



(58) 



4.4. 



The problem with the running Yukawa 
coupling 



We have been able to define a coupling constant at the 
scale M by direct calculation of the process f a + fb — > 
fa + fb with the scalar exchange in the t channel. It 
works as expected when dealing with this process near 
the renormalization point. Explicitly we see that 

-iM(f a + fb ^fa + fb)(q 2 ~ M 2 ) = 

(-i9l(M 2 ))-^{-igl(M 2 ))-i q l(M 2 ) 



. 3g a g b K 2 M 2 f J_ J_ 
8 16(4tt) 2 \\q 2 \ M 2 



■ 3g a gbK , f\q I 



Here we see that including the some of the loop effects in 
the coupling constant at the scale M 2 does a reasonable 
job of capturing the quantum effects in the neighborhood 
of this point. By construction, the residual quantum ef- 
fects vanish if we choose q 2 = —M 2 

However, now we consider the process f a + f a — > fb + fb 
with the scalar exchange in the s channel, with s = 
(pi + P2) 2 > 0, we see that that this definition of the 
running coupling has the wrong sign to correctly capture 
the effect of the quantum loops. Instead of incorpora- 
tion the majority of quantum effects at this scale gets 
the sign wrong so that the quantum effects are doubled. 
Explicitly we find 



-iM(f a + f a ^fb + fb)(s~M 2 ) 

i 
s 

„2 A/f2 /I 1 \ Q„ „. ^2 



(-igl(M 2 ))-(-igl(M 2 ))-i q l(M 2 ) 

3g a g b n 2 M 2 / 1 1 \ . -"V/,, '/,',/->- , / > 
16(4tt) 2 \s M 2 



.3g a gbK 



(59) 



(60) 

The process of absorbing these n 2 M 2 effects in the cou- 
pling constant has not described the kinematic depen- 
dence of the loops in this crossed reaction because the 
sign of q 2 has changed. The residual power-law effects 
are doubled at s = AI 2 , while the log effects associated 
with the running of 171 behave as expected. 

Alternatively we could have used the process f a + f a ^> 
fb + fb as our choice of renormalization procedure. In 
this case, the quantum correction and the running of the 
coupling g in the beta function Eq. [57] would have the 
opposite sign. This defintion would work fine for the s- 
channel process, but would then fail for the f-channel 
reaction. 

It is clear that we can have a sensible on-shell renormal- 
ization procedure which yields either sign for the running 
of the coupling, depending on whether we use the space- 
like or time-like reaction. Indeed there is a multiplicity of 
schemes that can generate a wide range of answers. It is 
also clear from this that no single scheme (either on-shell 
or off-shell) will be able to correctly categorize the lead- 
ing quantum effects in all processes, because the crossed 
reaction will have a result of the opposite sign. This is 
to be contrasted with the definition of running coupling 
in renormalizable theories, which is universally success- 
ful for all reactions independent of the renormalization 
scheme. 

The difference comes from the phenomenon of opera- 
tor mixing. In the case of renormalizable theories, the 
quantum effects are renormalizing the original operator. 
However in effective field theories the quantum effects are 
associated with an operator of higher dimension. We can 
by fiat take some of the higher order effect and build it 
into the original coupling - there are many ways to do 
this. However, we see that such a construction fails to 
account for the quantum effects in related processes in 
the correct way. 

There are further process dependence that is evident 
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in our results. If we were to consider a different on-shell 
reaction, say 4>-\-ij) — > </> + 4>, this would involve the basic 
(j>ipif! vertex with a different particle off-shell, in this case 
the fermion. Our result for the off-shell renormalization 
shows that this involves a different numerical factor, and 
equivalently the overall process involves different opera- 
tors in the Qi basis of Eq. 221 Again we see that because 
the loops renormalize higher order operators that do not 
lead to a universal running coupling for all processes. 



5. BRIEF COMMENTS ON THE LITERATURE 

The study of gravitational corrections to the running 
of gauge couplings was started by Robinson and Wilczek 
who calculated the one-loop contribution of graviton 
exchange to the f3 function of the Yang-Mills theory us- 
ing the background field method. Subsequently a series 
of authors used a variety of different methods to 

argue that the gravitational correction to the running 
of gauge couplings actually vanishes. Yet more recently, 
there have been further studies [l(| that again claim 
non-zero effects for the running of gauge couplings. Sim- 
ilar discrepancies are found in work that studies the non- 
gauge interactions of scalars and ferrnions 

EMU 

again 

using a variety of methods. 

The variety of answers found using different methods is 
itself a indication of the non-universality of gravitational 
corrections to running couplings. Because of the power 
counting for gravitational loop corrections, the true effect 
of gravitational loops is to renormalize a higher order 
operator and not the original vertex. Of course, in some 
schemes with a dimcnsionful cutoff, there is an additional 
renormalization of the original operator that is not seen 
when using dimensional regularization. However, this is 
a scheme dependent artifact and we comment on this 
below. 

In addition, most schemes in the literature do not con- 
sider the effects of high order operators nor of process 
dependence. By mimicking the effect of higher order op- 
erators through a redefinition of the original coupling, 
one can create the appearance of a running coupling. 
But unless this definition is universally valid, it is not 
a true running coupling. Our work shows the limitations 
of such definitions. Moreover, as we have shown, the 
operator basis plays a different role in off-shell methods 
versus on-shell methods. 



5.1. Comments on cutoff regularization 

We have performed our calculations using dimensional 
regularization. Real physics does not depend on the 
renormalization scheme, so that we could obtain equiva- 
lent results in any consistent regularization scheme that 
respects general covariance. However, because the grav- 
itational coupling is dimensionful, the cutoff can appear 



as a power in the renormalization procedure and this can 
create some confusion. 

Let us review how one would recover the results of di- 
mensional regularization in situations in which a cutoff is 
used. One imagines that the momentum range is divided 
into regions above and below an arbitrary cutoff A. The 
contribution of the region above the cutoff is of course 
unknown, and the effects are contained in the coupling 
constant <7 bare (A) which is parameterized by the separa- 
tion scale A. The theory below the cutoff is treated as 
a full field theory and the loops are cutoff by A. In the 
case of gravity this cutoff appears quadratically because 
of the dimensional gravitational coupling. The bare cou- 
pling and the loop correction are added together to yield 
the renormalized physical coupling 

gP hys. = ff barc( A ) + ^2^2 ^ 

which however is independent of the arbitrary cutoff A. 
It is the physical coupling which is equivalent to the cou- 
pling used in a dimensionally regularized scheme. In di- 
mensional regularization, the loop momentum runs over 
all energies and there is no need for a separation scale, 
nor for a divergence in the renormalization of the original 
operator. As expected then, physical effects are the same 
with either regularization scheme. 

Should one identify the dependence of the bare cou- 
pling <? bare (A) on the cutoff as the running of the coupling 
constant? The answer is clearly negative. That depen- 
dence cancels completely in physical observables. It also 
is independent of the kinematic dependence of scattering 
amplitudes so that it does not capture the real quan- 
tum effects of physical processes. This situation is differ- 
ent in the case of renormalizable field theories. In those 
cases the dependence on the cutoff is logarithmic, and by 
dimensional grounds the kinematic variables also enter 
into the logarithm, ln(A 2 /g 2 ). So even though the cutoff 
cancels in physical observables it nevertheless provides a 
guide to the kinematic dependence of the quantum cor- 
rections. The counterpart in dimensional regularization 
is the study of the 1/e and ln(fi 2 /q 2 ) terms which also 
mirror the kinematic dependence of loop effects. The fact 
that 1/e or ln(A 2 ) goes into the renormalization of the 
coupling also leads to the universality of these quantum 
effects, independent of the process being considered. 

We see that the quadratic cutoff dependence in gravi- 
tational corrections does not signal the appearance of a 
running coupling in physical processes but totally drops 
out as an unphysical artifact. This is why we focussed di- 
rectly on the kinematic dependence of physical processes 
and the influence of that on a coupling through the choice 
of renormalization procedure. This type of behavior is 
independent of the regularization scheme. 

6. DISCUSSION 

Effective field theory techniques expand a theory about 
zero energy, and the renormalization and higher momen- 
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turn dependence involves the operators that appear at 
higher order in the energy expansion. Here the content 
of the renormalization group is limited to the descrip- 
tion of the coefficients of the leading logarithmic loop 
corrections, but the couplings of the theory do not run 
with the energy scale. In contrast, there are many at- 
tempts to define running coupling constants in such the- 
ories by renormalizing the theory at a higher energy scale 
M . We have explored this process using off-shell and on- 
shell renormalization techniques, and have studied the 
subsequent utility of a running coupling in perturbativc 
scattering processes. This provides a well-defined setting 
for exploring the nature of these running couplings. 

It is of course possible to define running couplings in 
effective theories such as gravity. Indeed, there are a 
quasi-infinite number of ways of defining such couplings, 
so the question is not whether it can be done but rather 
whether it is useful and also to what extent it is univer- 
sal. The comparison standard is the definition of run- 
ning couplings in renormalizable theories where, despite 
minor scheme dependence, the results are highly useful 
and universal. 

Unfortunately, the corresponding definitions including 
gravity are in general seen not to be useful or univer- 
sal when describing physical reactions. The corrections 
go like n 2 q 2 , and in physical processes q 2 can be either 
positive or negative. For theories such as the Yukawa 
couplings, it is not possible to define a running coupling 
that is appropriate for both space-like and time-like pro- 
cesses. A theory such as <j) A is seen to be an exception. 
Gauge theories and gravity coupled to matter would be- 



have similar to the Yukawa theory. 

Instead, we see that in the perturbative realm the 
physics which would be described by a running cou- 
pling is accounted for by operator mixing. Higher di- 
mension operators with derivative coupling appropriately 
deal with both space-like and time-like processes. 

In our studies we also demonstrated how operators 
which vanish by field redefinition or by the equations 
of motion, such that they would normally be discarded 
from the operator basis, are nevertheless required when 
one performs the renormalization at an off-shell point. 
This is not surprising because the general theorems which 
demonstrate the irrelevance of such operators only ap- 
ply on-shell. For this reason the effects of these oper- 
ators continues to disappear when applied to physical 
processes. However, it points to the need of an expanded 
operator basis for field theory methods that involve off- 
shell renormalization. 
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